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We have investigated the imaginary part of the dielectric function Im(e) of the 
(113) 3x2 ADI reconstructed surface of sihcon. The calculations have been 
performed for a periodic slab within the plane-wave pseudopotential approach 
to the density-functional theory. The three diagonal components of Im(e) have 
been derived from the momentum matrix elements within the independent 
particle random phase approximation (IPRPA). In this article, the importance 
of the k-point convergence is figured out by inspecting k-point resolved spectra. 
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1. Introduction 

The Si(113) surface is one of the most stable high-index surfaces of Si ^. 
It is used as a substrate for the self-assembled growth of Ge nanowires 
and Ge as well as SiGe islands^"^. Since this surface is atomically smooth, 
ultrathin oxide films can be grown on Si(113)^'^, which is hence also dealt 
as a candidate for wafers and devices. 

The clean surface shows a 3x2 periodicity at room temperature 
which can transform into a 3x1 one by a phase transition induced 
by high temperatures or contamination'*'^. The most probable is the 
Si(113)3x2ADI (adatom-dimer-interstitial) reconstruction, which was pro- 
posed by Dabrowski et al.^° (see Fig. 1). Comparing ab-initio results for 
reasonable surface models, the ADI reconstruction has the lowest sur- 
face energy^^. Besides, an experimental confirmation is given by STM 
measurements^' However, still other models can not be ruled out^'^. 

A complementary study of the Si(113)3x2 surface can be done by the 
investigation of its optical properties, e.g., the reflectance anisotropy spectra 
(RAS). At the aim of performing a theoretical ab-initio investigation of the 
optical properties the first step is the calculation of the imaginary part of 
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Fig. 1. In the left panel, a top view {x-y plane) of the Si(113)3x2ADI reconstruction 
is shown. The right panel displays a top view of just two "double layers" of the bulk-like 
part of the slab. The Si atoms are dark (the interstitial Si atom light) and the hydrogen 
atoms light. 

the dielectric function Im(e). Results for Im(e) are presented in this work. 
2. Method 

The calculations have been performed within the plane-wave (PW) pseu- 
dopotential approach to the density- functional theory -'^^'-'^^ as implemented 
in ABINIT^^ and TOSCA^^. For the exchange-correlation energy the local- 
density approximation has been used^®'^^. The ground state convergence 
required 4 k Monkhorst-Pack^° points in the irreducible wedge of the Bril- 
louin zone (IBZ) as well as a kinetic-energy cutoff of 12 Ry (w 24000 PW). 
The slab used for the calculation contains 11 double layers (DL) of Si where 
the bottom surface is saturated with hydrogen atoms (see Figs. 1 and 4). 
The topmost 4 DL have been relaxed, and the remaining have been kept 
to bulk positions. 

The imaginary part of the dielectric function Im[e(ti;)] as a function of 
the energy uj has been derived by a summation of the matrixelements of the 
momentum operator^^. The sum has to be taken over the valence- and the 
conduction states (here: 175 and 185, respectively), and weighted k points, 
which in the present case can be taken just in one quarter of the Bril- 
louin Zone, which is the IBZ of the system. For the layer-by-layer analysis 
modified matrixelements have been used according to the prescriptions of 
reference 22 (see also references 23 and 24). Here, the effects of local fields, 
self energy, and excitons have been neglected. 



February 6, 2008 22:12 



WSPC - Proceedings Trim Size: 9in x 6in gaalnagy 



3 




1 2 3 1 2 3 

w (eV) w (eV) 



Fig. 2. Im[ea;{Lj)] (left) and Im[€y{u>)] (right) for various numbers of k points (assign- 
ment in the graph) for the low-energy range, where the surface states are found. The full 
energy range is shown in the insets. 

3. k point convergence 

For the calculation of Im[e(ijj)] (and consequently, for the calculation 
of all optical properties which are based on this) the convergence with re- 
spect to number of k points in the summation described above is essential, 
since the spectra are particularly sensitive to that. Working out of conver- 
gence can yield wrong difference spectra (e.g., RAS). In comparison, the 
convergence with respect to the number of bands is less demanding in the 
surface-related low-energy range. Thus, we performed a very careful check 
of the k-point convergence. 

Since the slab was chosen with the surface perpendicular to z, Im[e(ci;)] 
should be investigated for light polarized in x and y direction, since these 
functions arc the ingredients for difference spectra like the RAS. We show in 
Fig. 2 the Im[e2;(ti>)] and Im[ej^(aj)] for various sets of k points. For both po- 
larizations, at least 25 k-points are required. The convergence for Im[ey(w)] 
is faster than for Im[ea;(ijj)]. In particular, the differences between 4 and 25 
k points for the Im[ej:(a;)] is still significant. 

In order to improve the k-points grid in the most efficient way, we have 
performed a single-k analysis. This means, we have inspected the contribu- 
tions of each k point to the spectra separately, which are shown in Fig. 3. It 
is apparent that the changes for the Im[e2,(a;)] from one k to a neighbouring 
one, especially in x direction, are larger than for Im[ey(w)]. Thus, we have 
chosen 30 additional k points in between along the x direction. The spectra 
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Fig. 3. Single-k decomposition of the spectra using 25 k in the surface IBZ. Im[ex(aj)] 
is drawn with sohd and Im[ej^(ii;)] with dashed hnes. The coordinates of the k points are 
denoted in the inset. The spectra are organized as the k points in the surface IBZ, where 
the r point is close to the low-left corner of the figure. 



based on the resulting 55 k points (not shown in Fig. 2) demonstrates, that 
convergence already had been achieved with 25 k points. This strategy al- 
lows one to avoid a "blind" increase of the k-point grid, which can be very 
time-consuming. 



4. Anisotropy at T 

In Fig. 2 a particular feature is visible: in the contribution coming from the 
F point, a strong anisotropy at the medium-high energy range, between the 
X and the y polarization, occurs. This anisotropy at F appears in the energy 
range of bulk-bulk transitions. In fact, performing a layer-by-layer analysis 
the spatial localization of such anisotropic contributions is clearly found to 
be in the bulk (internal) DLs, as shown in Fig. 4, together with a side view 
of the slab cell where the DLs used for each spectra are marked. We have 
chosen always pairs of DLs, e.g., L07L08 corresponds to the spectra with 
contributions from DL 7 and DL 8, where the DLs are counted from the top 
to the bottom. Of course, bulk silicon is isotropic, which results also from 
calculations using a sufficient amount of k points. The bulk anisotropy at 
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Fig. 4. Layer-by-layer decomposition of the spectra and complete spectra containing 
all layers (full) for the T point calculation (right). The spectra correspond to the layers 
denoted in the inset and marked at the structure (left). Im[ej;(aj)] is drawn with solid 
and Im[ej,(u)] with dashed lines. In addition, the full-slab spectra are drawn with dotted 
lines. DL LOO corresponds to the vacuum at the top of the slab and DL12 to the hydrogen 
at the bottom. 

the r point here is just due to the use of a low-symmetry supercell. 

Fig. 4 shows that this anisotropy is present at each pair of DLs in the 
bulk region of the slab, where the spectra for bulk region, i.e. L09L10, 
L07L08, and L05L06, are nearly the same. The explanation can be found 
by inspecting the geometry of one pair of DLs only, which is shown on the 
right in Fig. 1. An highly anisotropic "chain structure" oriented along the 
X axis can be seen. This is the reason for the strong signal appearing at F 
(the wavefunctions are in phase and the wavelength of light is much larger 
than the cell dimension). 

5. Summary 

In summary, we have presented ab-initio results for the imaginary part of 
the dielectric function for the Si(113)3x2ADI surface. At this example we 
have shown the possibility of accelerating the convergence tests with respect 
to k points by performing a single-k analysis. We have applied a layer- 
by-layer analysis allowing us to explain a strong anisotropic contribution 
coming from the F point. The results described here are useful for further 
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calculations, e.g., the RAS spectra, which can be obtained from Im[e(a;)]. 
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